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Abstract 



The dynamics of a coupled Bose-Einstein condensate involving trapped 
atoms in two quantum states is studied using the time-dependent Gross- 
Pitaevskii equation including an interaction which can transform atoms from 
one state to the other. We find interesting oscillation of the number of atoms 
in each of the states. For all repulsive interactions, stable condensates are 
formed. When some of the atomic interactions are attractive, the possibil- 
ity of collapse is studied by including an absorptive contact interaction and 
a quartic three-body recombination term. One or both components of the 
condensate may undergo collapse when one or more of the nonlinear terms 
are attractive in nature. 

PACS Number (s): 03.75.Fi 
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I. INTRODUCTION 



Recent experimental success |T|-|5| in Bose- Einstein condensation (BEC) at ultralow tem- 
perature in dilute trapped alkali metal and hydrogen atoms has intensified theoretical studies 
of the inhomogeneous, weakly interacting dilute Bose gas and its condensation ||6|-[r7|. At 
zero temperature such a system is supposed to be fully condensed and described by the 
mean-field Gross-Pitaevskii (GP) equation [|g|TH],|Tg. 

Two most interesting phenomena in BEC are the collapse for attractive interaction and 
the formation of coupled atomic condensate. For attractive atomic interaction the 
condensate is stable for a maximum critical number of atoms. In the presence of an external 
source of atoms when the number increases beyond the critical number, due to interatomic 
attraction the radius of BEC tends to zero and the central density tends to infinity. Con- 
sequently, the condensate collapses emitting atoms until the number of atoms is reduced 
below the critical number and a stable configuration is reached. Then the condensate can 
grow again and a series of collapses can take place, which was observed experimentally in 
the BEC of ^Li with attractive interaction f^-^]. 

There has also been experimental realization of BEC involving atoms in two different 
quantum states p0|-p2|. In one experiment ®^Rb atoms formed in the \F = 1, m = —1) and 
|F = 2, m = 1) states by the use of a laser served as two different species, where F and m 
are the total angular momentum and its projection [^^. In another experiment a coupled 
BEC was formed with the ^''Rb atoms in the \F = l,m = —1) and |F = 2,m = 2) states 



1^ , |22[| . It is possible to use the same magnetic trap to confine atoms in two quantum states 



and this makes these experiments technically simpler compared to a realization of BEC with 
two different types of atoms requiring two independent trapping mechanisms. 

Here we study theoretically the dynamics of a coupled BEC composed of two quantum 
states 1 and 2 of an atom using the coupled time-dependent GP equation. We include an 
interaction term in the Hamiltonian that allows for a coherent boson-number-conserving 
nondissipative transition of atoms from one quantum state to another. This interaction 
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simulates the action of a laser which can been used experimentally to transform one species 
of atoms to another. This model exhibits rich and novel phenomenology which we study 
using the numerical solution of the coupled GP equation, which can generate oscillation in 
the number of atoms in a specific state in the presence of the interaction that transforms 
one spices of atoms to another. Also, depending on the nature of the atomic interactions 
one can have three possibilities of BEG: (i) stable condensate of both components for all 
repulsive interactions, (ii) collapse of both components for all attractive interactions, and 
(iii) a stable condensate for one component and collapse in the other. 

The collapse for attractive interaction(s) in a coupled condensate is studied by intro- 
ducing absorptive contact interactions responsible for a growth in the atomic numbers from 
external source in addition to imaginary three-body quartic interactions leading to recom- 
bination loss In the presence of these imaginary interactions the GP equation does 
not conserve the overall number of atoms. If the strengths (and signs) of these imaginary 
interactions are properly chosen, one can have collapse in one or both components of the 
condensate. The solution of the GP equation could produce a growth of the condensate 
with time for atom numbers less than the critical numbers for collapse. Once they increase 
past the critical numbers, the three-body interaction takes control and the numbers sud- 
denly drop below the critical level by recombination loss signaling a collapse Then 
the absorptive term takes over and the individual numbers start to increase. This continues 
indefinitely showing an infinite sequence of collapses for one or both components. 

We motivate this study by considering two possible atomic states 1 and 2 of ''Li (attractive 
interaction) and of ^^Rb (repulsive interaction) whenever possible. In the case of ''Li the 
interaction in state 1 is taken to be attractive which is responsible for collapse. We study 
the collapse with different possibilities of attraction and repulsion between atoms in states 
1 and 2 of ^Li. 
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II. COUPLED GROSS-PITAEVSKII EQUATION 



We consider the following Hamiltonian to describe the coupled BEC allowing for the 
possibility of transforming atoms from one state to another 



H = Hq + Vi + V2, 



(1) 



where the kinetic and potential energies are 



2m 



(2) 



(3) 



V2 = xy dr[ijl{r)ilJ2ir) + ^lir)Mr}], (4) 

where ipj{r) and ^j{r) are the field operators for annihilation and creation of the bosonic 
atom of mass m in state j = 1,2, V^^\r) is external trapping potential, Vi is the usual 
nonlinear potential between the atoms interacting via a contact interaction of strength gij = 
Anh'^aij/m between two atoms in states I and j, aij is the scattering length of atoms in states 
/ and j, and V2 is the potential which allows for the laser-induced transition of an atom from 
state 1 to 2 and vice versa via a contact interaction of strength x- 

Coupled mean- field GP equations are obtained by replacing the field operators ipj (^) by 
the wave functions ipj{r), which are the expectation values of the field operators, in the 
Heisenberg equation [O]: 



(5) 



where the time (r) dependence is exphcitly shown and the result is p3|-p5| 
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V',(r,r) + 5]x(l-5,7)^z(r,r) = 0, 
1=1 



(6) 



where 6ji is the Kronecker delta and rii is the number of condensed atoms in state /. Here 
we have speciahzed to the radially symmetric case and have used the harmonic oscillator 
trapping potential V^^{r) = cjmuj'^r'^ /2 where uj is the trap frequency and Cj is a parameter 
introduced to vary independently the trap potential in each quantum state. The normaliza- 
tion condition for x = is given by 
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Att / r'dr\iljJr,T)\' = 1. (7) 

JO 

As in Refs. it is convenient to use dimensionless variables defined by x = y^r/a-^^ 



, and t = Tu, where aj^^ = ^Jh/{miu), and (j)j{x,t) = (r, r)(^/27ra^^)^/^ . In terms of 
these variables Eq. (|) becomes |]rT| 



dx' 4 x2 x^ 
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(Pjix, t) + J2 ^(1 - t) = 0, (8) 

1=1 

where rj = x/{^^) and riji = 2\/2niaji/a-^Q could be negative (positive) when the corre- 
sponding interaction is attractive (repulsive). In Eqs. (§) we have introduced a diagonal 
absorptive i'jj and a quartic three-body term —iC^j\(j)j{x, t)\'^/x'^ appropriate to study collapse 
0J^. A nonzero allows the possibility of the absorption of atoms from an external source 
into the condensate, whereas a nonzero C,j allows for the possibility of ejection of atoms 
from the condensate due to three-body recombination. The total number of atoms in the 
condensate is not conserved for nonzero 7-,- and ^j. 

The fluctuation of the number of atoms in the two states of the condensate is best studied 
via the quantities 
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N,= / \<j),{x,t)\'dx, J = 1,2. (9) 
Jo 

In the presence of a general absorptive interaction rj ^ 0,7j 7^ 0, and 7^ 0, the reduced 
number for the two components of the condensate are given by riuNi and n22N2, the quan- 
tities A^^i and N2 carry the information about time evolution of the number of the two 
components. The actual number of atoms in the two components are njNj, j = 1, 2. 



For 7j = = ''7 = 0, the number of atoms in the two components are conserved 
separately. This condition leads to Nj = 1. For 7^- = = and 1] 0, the numbers of 
atoms in the two components are not conserved separately Nj 7^ 1; but the total number 
is conserved. In the presence of absorptive interaction 7^ 7^ and/or 7^ 0, there is no 
conservation of even the total number of atoms. 



III. NUMERICAL RESULTS 



A. Computational Procedure 



To solve Eqs. numerically the proper boundary conditions as x ^ and 00 are 
needed. For a confined condensate the asymptotic form of the physically acceptable solution 
is given by \imx^oo\4>j{x,t)\ ~ exp(— At x = one has the regularity condition 
0,(O,t)=O. 

Next we discretize Eqs. (|) in both space and time by using a space step h and time step 
A with a finite difference scheme using the unknown (0^)^ which are approximation to the 
exact solution (j)j{xp,tk) where Xp = ph and = kA. After discretization we obtain a set 
of algebraic equations which could then be solved by using the known asymptotic boundary 
conditions |]6|,|7|,p!Tl| . The time derivative in these equations involves the wave functions at 
times tk and t/c+i- To discretize Eqs. (§) we express the wave functions and their derivatives 
in terms of the average over times tk and t^+i and the result is the following Crank- 
Nicholson-type scheme 
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1=1 



where j = 1,2. 



Considering that the wave function components ipj are known at time tk-, Eqs. ( [lOD 
involve the unknowns — {4>j)p^^ and (0j)pli at time tfc+i- In a lattice of points 

Eqs. ([T0| ) represent a tridiagonal set for p = 2, 3, (A^ — 1) for a specific component This 
set has a unique solution if the wave functions at the two end points {4>j)i^^ and are 
known . These values at the end points are provided by the known asymptotic boundary 
conditions. 

To solve Eqs. (|l^) we employ space step h = 0.0001 with xmax < 15 and time step A 
=0.05. After some experimentation we find that these values of the steps give good conver- 
gence. The iteration is started with the known normalized (harmonic oscillator) solution of 
Eqs. obtained with riji = 7^- = = rj = 0. The nonlinear parameters nji are increased 
by equal amounts over 1000 time iterations starting from zero at t = —50 until the desired 
final value is reached at t = 0. During the iteration we keep C,j = 7i = and attribute a 
small numerical value to rj (< 0.2). The resulting solution is the ground state of the con- 
densate corresponding to the specific nonlinearity and built up with this particular value of 
1]. Then to study the dynamics we allow the system to evolve by continuing the iteration 
with A = 0.05 and h = 0.0001 for positive t, but maintaining rj and Uij at the constant 
final values attained at t = 0. To study the dynamics for positive t for attractive atomic 
interaction, the absorptive contact (7^) and the three-body interactions (^j) are switched on 
at t = 0. This will clearly show the possible collapse of the system. 



B. Stable Condensate: Repulsive Interaction 

First we consider the stationary solution of Eqs. (§) with 7^- = C,j = in the repulsive 
case. We now choose the following two sets of of parameters before the actual calculation: 

= = 10, ni2 = = 5, Ci = 1, C2 = 0.25, {0)1] = 0.075 and (b) t] = 0.1. This will 
show the effect of the feedback between the two atomic states. In this case all interactions 
are repulsive corresponding to positive sign of all Uji = 2^/2njaji/a-^^. 

Although the parameters above are in dimensionless units, it is easy to associate them 
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with a physical problem of experimental interest in the case rj = 0. For example, for the 
mixture of |F = l,m = —1) and |F = 2,m = 1) states of '^^Rb, the ratio of scattering 



lengths a|i,-i)/ct|2,i) = 1-062 ||2T|. If we label state |1, —1) by 1 and |2, 1) by 2, and consider 
aii/flho ~ ~ 0.002, then riu = ^22 = 10 corresponds to rii ~ ^2 — 1770 for r] = 0. 

Hence these models can simulate the actual experimental situation composed of two states 
of 87Rb. 

In Fig. 1 we plot the wave functions \(f)j{x,0)\/x for the two components calculated 
with two nonzero values of rj above. In Figs. 2 (a) and (b) we plot the quantities A^i 
and N2 for these models. When rj = 0, Ni = N2 = 1. As a nonzero value of rj is taken, 
the quantities Nj oscillate with time as can be seen from Figs. 2 (a) and (b). When A^i 
increases, N2 decreases and this denotes a net transformation of atoms from state 2 to 1 and 
vice versa. These transformations lead to interesting oscillation in Nj and consequently, in 
the actual number of atoms in the two states. For a small rj these oscillations are small and 
as rj increases these oscillations increase in amplitude. These oscillations clearly show the 
continued dynamical conversion of one species of atoms to another in the otherwise stable 
condensate. 



C. Collapse: Attractive Interaction 

Now we study the simplest case of collapse by taking only the interaction between the 
atoms in state 1 to be attractive corresponding to a negative an. All other scattering 
lengths — a22 and (= 021) — are taken to be positive. Quite expectedly, here the first 
component of the condensate could experience collapse. We also consider the collapse in 
both components when all scattering lengths are negative. 

The calculation is performed in these cases with the following two sets of parameters (a) 
rill = —3.0346, ^22 = 4, = n2i = 1, Ci = 0.25, C2 = 4, r/ = 0.2 and (b) riu = —0.67, 7222 = 
— 1.4, ni2 = — 0.42,'n,2i = —0.41, ci = 4, C2 = 0.25, rj = 0.1. The wave function components 
at t = in these cases are plotted in Fig. 3. 
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The above parameters are in dimensionless units and one can associate them with a 
problem of experimental interest for r/ = in terms of two states of ^Li. In case of model 
(a) above we consider the state 1 to be the ground state of ''Li with attractive interaction as 
in the actual collapse experiment with |ctii|/%Q — 0.0005 As riu = 2v^r2i|aii|/aj^Q, 

this corresponds to a boson number ni ~ 2145. This number is larger than the maximum 
number of atoms permitted in the BEC of a single component ^Li which is about 1400 . 
The state 2 could be one of the excited states of ^Li with repulsive atomic interaction. The 
presence of the second component with repulsive interaction allows for the formation of a 
stable BEC with more ^Li atoms in quantum state 1 than allowed in the single-component 
BEC. We find from Fig. 3 that (pi (full line) is very much centrally peaked compared to 
02 (dashed line) in this case. This corresponds to a small radius and large central density 
for the first component denoting an approximation to collapse. If the number Ui is slightly 
increased beyond 2200 the first component of the condensate wave function becomes singular 
at the origin and could experience collapse. 

In case of model (b) we take the atomic interaction in both states of ^Li to be attrac- 
tive corresponding to negative scattering lengths. For |aii|/aj^Q ~ 0.0005 as in the actual 
experiment one has ni = 470 and n2 — 990 for rj = 0. The total number of particles 
in this case is roughly 1460, which is close to the critical number 1400 found in the actual 
experiment of collapse in ''Li. Both wave-function components could become singular in this 
case as all possible interactions are attractive. The small extension of the wave-function 
components in space in this case as in Fig. 3 denote an approximation to collapse. 

Although the collapse of the coupled condensates could be inferred from the shape of 
the stationary wave functions of Fig. 3 (sharply peaked centrally with small radii), we also 
study the dynamics of collapse from a time evolution of the full CP equation (|]) in the 
presence of an absorption and three-body recombination, e.g., for 7^- 7^ and C,j 7^ as in 
the uncoupled case 0. The general nature of time evolution is independent of the actual 
values of jj and employed provided that a small value for and a relatively larger one 
for 7j are chosen p. 
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First we consider model (a). For a dynamic description of this problem we take for t > 
7i = 72 = 0.04, ,^1 = 0.01 and ,^2 = 0.002 and allow the solution of the GP equation to 
evolve in time using Eq. (0). In the presence of these absorptive interactions the total 
number of atoms is not a constant of motion. The time evolution of A^^i and N2 is shown 
in Fig. 4 (a). There is small oscillation in Nj for t < signifying a transfer of one type of 
atoms to the other. For t > these oscillations are almost undetectable. Because of the 
attractive interaction, the component 1 could experience a succession of collapse. This is 
manifested in the successive growth and decay of A^^i. The condensate corresponding to the 
second component controlled by the repulsive interaction keeps on growing, manifested by 
the continuous growth of N2 with time. 

Next we consider the time evolution of the GP equation for model (b) with the following 
absorptive parameters: 71 = 0.04,72 = 0.03, ^1 = ^2 = 0.005. In this case as all interactions 
are attractive, both components can experience collapse. In Fig 4 (b) we plot the evolution of 
A^^i and N2- The oscillation of A^^i and N2 is clearly visible for t < 0. For t > 0, the absorptive 
interactions take control and the total number of atoms start to increase. Consequently, the 
dynamics of collapse is more visible than the oscillation in A^^i and N2- Both A^i and N2 is 
found to exhibit a succession of growth and decay corresponding to collapse. Hence we find 
that one can have a dynamical collapse of one or both components of the BEG in the case 
of a coupled condensate. 

IV. CONCLUSION 

To conclude, we studied the collapse in a trapped BEG of atoms in states 1 and 2 
using the GP equation when some of the atomic interactions are attractive. We include 
an interaction in the Hamiltonian which allow for a spontaneous transformation of atoms 
from state 1 to 2 and vice versa. Experimentally, this can be achieved by a laser and makes 
the present study richer. In case of all repulsive atomic interactions one can have a stable 
condensate of both components. When some of the atomic interactions are attractive one 



10 



or both components of the condensate could experience collapse. The collapse could be 
predicted from a stationary solution of the GP equation. The time evolution of collapse is 
studied via the time-dependent GP equation with absorption and three-body recombination. 
The number of particles of the component (s) of BEG experiencing collapse alternately grows 
and decays with time. With the possibility of a detailed study of a coupled BEG, the results 
of this study could be verified experimentally in the future. 

The work is supported in part by the Gonselho Nacional de Desenvolvimento Gientifico 
e Tecnologico and Fundagao de Amparo a Pesquisa do Estado de Sao Paulo of Brazil. 
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Figure Caption: 

1. Wave function components \(l)i{x,0)\/x (full line) and |02(a^, 0)|/a: (dashed line) vs. x 
for two coupled GP equations with nu = ^22 = 10, ni2 — n2i — 5, ci = 1, C2 = 0.25, {a)r] — 
0.075 and (b) 77 = 0.1. 

2. A^i (full hne) and N2 (dashed hue) vs. t for models (a) and (b) of Fig. 1. 

3. Wave function components \4>i{x,0)\/x (full line) and \4>2{x,0)\/x (dashed line) vs. x 
for two coupled GP equations with (a) nu — —3.0346,7122 = 4, ni2 — n2i — 1, ci — 0.25, 
C2 = 4, 77 = 0.2 and (b) nn = -0.67, 7122 = -1-4, 7ii2 = -0.42, 7121 = -0.41, ci = 4, 
C2 = 0.25, 7] = 0.1. 

4. A^i (full hne) and N2 (dashed hne) vs. t for models (a) and (b) of Fig. 3. 
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Figure 1 




Figure 2 (b) 
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